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ABSTRACT 

A simple, physically consistent model has been proposed that seeks to explain in a 
unified way the X-ray spectra and rapid X-ray variability of the so-called Z sources and 
other accreting neutron stars in low-mass systems. Here we summarize the results of 
detailed numerical calculations of the X-ray spectra of the Z sources predicted by this 
model. Our computations show that in the Z sources, photons are produced primarily by 
electron cyclotron emission in the neutron star magnetosphere. Comptonization of these 
photons by the hot central corona and radial inflow produces X-ray spectra, color-color 
tracks, and countrate variations like those observed in the Z sources. 
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1. INTRODUCTION 

The six known Z sources (Hasinger & van der Klis 
1989) are among the most luminous low-mass X-ray 
binaries. They are unique in tracing out Z-shaped 
tracks in X-ray color-color diagrams on time scales of 
minutes to hours (see van der Klis 1989; van der Klis 
& Lamb 1995). The tracks are thought to be pro- 
duced by variations in the mass accretion rate, which 
appears to increase monotonically as a Z source moves 
from its horizontal to its normal and then to its flaring 
branch (see van der Klis & Lamb 1995). Power spec- 
tra of Z-source brightness variations show two types 
of quasi-periodic oscillations (QPOs) and three types 
of aperiodic flickering ("noise"). The properties of 
the QPOs and flickering vary systematically with the 
position of a source on its Z track. 

The X-ray spectra of the Z sources and the origin of 
their Z-shaped tracks has been a puzzle. As the preci- 
sion of X-ray spectral measurements has improved, a 
variety of increasingly complicated mathematical ex- 
pressions and physical models has been used in at- 
tempts to obtain statistically acceptable fits to the 
spectra of Z sources and other accreting neutron stars 
in low-mass systems (see, e.g., Swank & Serlemitsos 
1985; Hirano et al. 1984; Mitsuda et al. 1984, 1989; 
Czerny, Czerny, & Grindlay 1986; White et al. 1986; 
Vacca et al. 1987; White, Stella, & Parmar 1988; Pon- 
man, Foster, & Ross 1990; Christian & Swank 1995). 
No attempt was made in these works to relate the 
spectral properties of a source to its rapid X-ray vari- 
ability, and many of the models are not physically 
consistent (see Lamb 1986; White et al. 1988). In 
these early works, models were fit to X-ray spectra 
averaged over hours, and the origin of the Z tracks 
of the Z sources was not addressed (in some cases the 
existence of the Z tracks was unknown at the time the 
work was done) . 

Several authors have shown that the tracks of one 
or two Z sources can be reproduced by varying the pa- 
rameters in various complicated expressions (Schulz, 
Hasinger, & Triimpcr 1989; Hasinger et al. 1990; 
Hoshi & Mitsuda 1991; Honma & Mineshige 1993), 
but the required variations of the parameters are typ- 
ically erratic and not motivated by physical argu- 
ments. Schulz & Wijcrs (1993) showed that suitable 
variations of the parameters in a model which assumes 
that blackbody emission from the neutron star sur- 
face is Comptonizcd by electrons in an extensive hot 
corona can reproduce the observed color-color tracks 



of four of the Z sources. However, the model is not 
physically self-consistent and the inferred variations 
of the parameters along the Z tracks are erratic and 
unrelated to the other properties of the source. 

A model that seeks to provide a unified, physically 
consistent explanation of the X-ray spectra and rapid 
X-ray variability of the Z sources as well as other 
low-mass binary systems containing accreting neutron 
stars has been proposed by Lamb (1989, 1995). In 
this 'unified model' the Z sources are rapidly spin- 
ning neutron stars with magnetic fields ~ 10 8 -10 9 G, 
accreting matter from their companions at ^0.5-1.1 
times the Eddington critical rate Me via a Keplerian 
disk and an approximately radial inflow that origi- 
nates in a corona above the inner disk. The radial 
flow supplies ~ 20% of the total mass flux. Interac- 
tion of the accretion disk and radial flow with the 
^3 x 10 6 cm neutron star magnctosphcrc produces 
a hot central corona (HCC) around the magneto- 
sphere. Soft (~ 0.5-0.8 keV) photons produced within 
the magnetosphere by high-harmonic cyclotron emis- 
sion and other processes are Comptonized by elec- 
trons in the HCC and the radial inflow. The interac- 
tion of radiation with the radial flow keeps its tem- 
perature close to the local Compton temperature (~ 1 
keV). The result is a relatively flat X-ray spectrum 
that cuts off near the ~ 10 keV temperature of the 
electrons in the HCC. Preliminary results (see Lamb 
1989) showed that the X-ray spectra and colors pre- 
dicted by the unified model are roughly consistent 
with the observed spectra and colors of the Z sources. 

When the accretion rate is ~ 0.5-0.9 Me, interac- 
tion of the small magnetosphere with the accretion 
disk produces a quasi-periodic luminosity oscillation 
at the beat frequency between the stellar spin fre- 
quency and the orbital frequency of clumps near the 
inner edge of the disk (Alpar & Shaham 1985; Lamb 
et al. 1985; Shibazaki & Lamb 1987). This is the 
so-called horizontal branch oscillation (HBO). When 
the mass accretion rate rises above ~ 0.8 Me, a soft, 
global, radiation-hydrodynamic mode of the radial in- 
flow becomes weakly damped and is therefore excited 
by fluctuations in the accretion flow, causing the X- 
ray spectrum to oscillate quasi-periodically at a fre- 
quency approximately equal to the inverse of the in- 
flow time from the outer radius of the flow (Lamb 
1989, 1995; Fortncr et al. 1989, 1995; Miller & 
Lamb 1992). At still higher accretion rates, nonra- 
dial modes develop (Lamb 1989, 1995; Miller & Park 
1995). This is the so-called normal/flaring branch os- 
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dilation (N/FBO). 

In this Letter we summarize detailed numerical cal- 
culations of the X-ray spectra of the Z sources based 
on the unified model. Our results provide a simple, 
physically consistent explanation of the spectra of the 
Z sources, including their Z tracks, as well as other 
previously unexplained properties of these sources. A 
more detailed account of the spectral model and fur- 
ther applications to the Z sources and other types of 
accreting neutron stars in low-mass systems will be 
presented elsewhere. 

2. SPECTRAL MODEL 

We have computed the spectrum of the radiation 
produced in the magnetosphcrc, taking into account 
resonant and nonresonant scattering as well as cy- 
clotron absorption and emission (Hartmann ct al. 

1988) and the effects of induced processes, which typi- 
cally dominate spontaneous processes by factors ~ 20. 
We assume that the intrinsic stellar magnetic field 
is dipolar and compute the variation of the mag- 
netic field with radius, treating the magnetospheric 
boundary as spherically symmetric. The spectrum of 
the radiation that emerges from the magnetosphere 
depends strongly on the field strength and electron 
temperature and more weakly on the electron den- 
sity. For stellar magnetic fields ~ 10 9 G, electron tem- 
peratures ~ 5—15 keV, and electron column densities 
~ 10 19 -10 20 cm~ 2 , the spectrum is roughly triangu- 
lar, following the Planck curve until somewhat below 
E max ~ 0.1-1 keV and then falling steeply as the 
thermalization length becomes larger than the mag- 
netosphere. 

Photons emerging from the magnetosphere into the 
HCC are Comptonized by the hot electrons there, los- 
ing energy to electron recoil and gaining energy from 
electron thermal and bulk motion (free-free absorp- 
tion and emission are negligible at the densities, tem- 
peratures, and photon energies of interest). The qual- 
itative shapes of Z-source X-ray spectra (Schulz et al. 

1989) suggest that T" cc , the electron temperature in 
the HCC, varies from 5 to 15 keV, depending on the 
source and its position on its Z track. The shapes of 
Z-source X-ray spectra and the low upper limits on 
the amplitudes of any periodic X-ray brightness oscil- 
lations (Vaughan et al. 1994) require that t" s cc , the 
electron scattering optical depth of the HCC, be ~ 5-7 
(Brainerd & Lamb 1987; Lamb 1989, 1995). Because 
t" s cc is tightly constrained by these considerations, 



we set it equal to 6 in the present calculations. 

Since T" cc is much greater than E max , the spec- 
trum of the radiation that emerges from the HCC 
above 1 keV is a slightly rounded power law that 
cuts off at ~T e HCC (see Rybicki & Lightman 1979). 
The exponent of the power law depends primarily 
on y HCC = (4k B T e HCC /m e c 2 ){T? s cc ) 2 , the Compton y 
parameter in the HCC. The mean energy of the pho- 
tons leaving the HCC is a factor ~ e v greater than 
the mean energy of the photons entering the HCC, so 
the luminosity of the HCC is L w e v L m , where L m is 
the luminosity of the magnetosphere. 

Interaction of the radiation from the HCC with 
the radial flow does not affect the luminosity (the 
thermal content of the flow is negligible and hence 
the temperature of the electrons in the radial flow 
is equal to the local Compton temperature; see Lamb 
1989, 1995), but it does affect slightly the shape of the 
spectrum, increasing the curvature of the spectrum at 
~ 2-5 keV. The qualitative shapes of Z-source X-ray 
spectra constrain r^f , the electron scattering optical 
depth of the radial flow, to be in the range 2-10; the 
radiation-hydrodynamic model of the N/FBO further 
constrains t*/ to be ~ 6 near the middle of the nor- 
mal branch (Fortner et al. 1989, 1995; Miller & Lamb 
1992). 

In computing the X-ray spectra presented here, we 
approximated the spectrum of the radiation coming 
from the magnetosphere by a triangular function that 
mimicks the results of our detailed numerical compu- 
tations and calculated the effects of Compton scatter- 
ing in the HCC and radial flow using the approach of 
Miller & Lamb (1992), which treats these regions as 
uniform. The finite size and spherical geometry of the 
HCC and radial flow are taken into acount by specify- 
ing appropriate distributions of photon escape times. 
The effect of the converging flow was neglected. For 
simplicity we assumed that the electrons in the mag- 
netosphere and HCC have the same temperature. To 
compare our calculated spectra with spectra observed 
with EXOSAT, we used the interstellar absorption 
cross-section of Morrison & McCammon (1983) and 
the appropriate energy response matrix of the EX- 
OSAT ME detector (see Turner, Smith, & Zimmer- 
man 1981). 

3. RESULTS 

The fundamental spectral parameters of the uni- 
fied model are the mass M, radius R, and magnetic 
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field B of the neutron star, and the mass accretion 
rate M. The X-ray spectrum observed at Earth also 
depends on two additional 'external' parameters: the 
distance D and the neutral hydrogen column density 
Nh between the Earth and the source. Of these six 
parameters, only M is expected to vary significantly 
on the time scales of interest here. 

In practice, some of the physical processes that are 
important near the neutron star are fairly well under- 
stood whereas others are not. For this reason we work 
with the three secondary quantities that directly enter 
our computation of the X-ray spectrum observed by 
EXOSAT, namely, the peak energy -E max of the spec- 
trum emerging from the magnetospherc, the Comp- 
ton parameter y HCC of the HCC, and the scattering 
optical depth r*f of the radial flow. As discussed 
in § 2, these parameters are tightly constrained by 
the qualitative shape of Z-source X-ray spectra and 
the rapid X-ray variability of these sources. Once the 
three secondary spectral parameters are determined, 
approximate physical relations can be used to esti- 
mate the corresponding values of B and M, thereby 
checking the consistency of the unified model. We de- 
termine the external constant N H from observations 
other than the ones with which we compare our pre- 
dicted spectra. 

We specify the position of a Z source on its track 
by a rank number, defined to be 1.0 at the junction 
of the horizontal and normal branches and 2.0 at the 
junction of the normal and flaring branches (see van 
der Klis & Lamb 1995). Figure 1 shows how E max , 
y HCC , and T^f are expected to vary with rank num- 
ber as M increases from 0.7 Me to 1.1 Me, m the 
unified model (see below). Figure 2 compares the re- 
sulting color-color and hardness-countrate tracks with 
Cyg X-2 tracks and the predicted countrate spectrum 
at rank number 1.5 with the countrate spectrum ob- 
served near the middle of the normal branch. The 
agreement is quite good, given the simplicity of the 
model. We now briefly describe the physical reasons 
for the changes in the spectral parameters shown in 
Figure 1, and explain why these changes produce the 
Z tracks shown in Figure 2. 

Horizontal branch. — As M increases from ~ 0.7 Me 
to ~0.9Mb, the magnetosphere is compressed, but 
£ max remains nearly constant because it is most sensi- 
tive to the magnetic field at ~ 2-3 stellar radii, which 
remains almost unchanged. The cyclotron luminosity 
increases more slowly than the accretion luminosity, 
so the electron temperature rises in order to main- 



tain thermal balance. The resulting increase in y HCC 
tends to increase both the soft and hard colors. The 
increase of the mass flux in the radial flow causes t* s f 
to increase, which tends to increase the countrate in 
the 2-5 keV energy range and hence to increase the 
soft color and decrease the hard color. The net ef- 
fect is to produce an approximately horizontal track 
in color-color and hardness-countrate diagrams. 

Normal branch. — As M increases from ~ 0.9 Me to 
~ 1.0 Me, the magnetosphere is further compressed. 
The increase in the mass flux and the rapid reduction 
of the inflow velocity caused by the increase in the 
outward radiation force cause n e to increase by a fac- 
tor ~ 8, causing E max to increase by ~ 10%. As a re- 
sult, the cyclotron luminosity increases faster than the 
accretion luminosity, and hence the electron temper- 
ature must fall in order to maintain thermal balance. 
The resulting decrease in y HCC causes the spectrum 
to become steeper, decreasing the soft and the hard 
colors. Meanwhile, t" s f continues to rise. As y HCC de- 
creases, fewer cyclotron photons are scattered up into 
the EXOSAT energy range and hence the countrate 
decreases, even though the luminosity is increasing. 
The result is a diagonal track running from upper 
left to lower right in both color-color and hardness- 
countrate diagrams. 

Flaring branch. — The character of the accretion 
flow when M exceeds Me is uncertain, but is proba- 
bly aspherical and time-dependent, with matter flow- 
ing inward in some regions and outward in others 
(Lamb 1989, 1995). In such a situation, radiation 
escapes primarily through the regions of lower opti- 
cal depth and the magnetosphere expands as M in- 
creases, decreasing _E max and increasing y HCC ; t" s f 
continues to increase. The result is a flaring branch 
that is nearly parallel to the normal branch but ex- 
tends to a slightly greater soft color. 

4. DISCUSSION 

Our results show that the unified model, which 
predicts that self-absorbed cyclotron emission in the 
neutron star magnetosphere is Comptonized by elec- 
trons in a hot central corona and radial inflow, pro- 
duces X-ray spectra that agree with those of Cyg X-2. 
Variation of M from 0.7 Me to 1.1 Me produces Z- 
shaped tracks in X-ray color-color diagrams (CDs) 
and hardness-countrate diagrams (HCDs) that agree 
with the tracks made by Cyg X-2. 

Analysis of archival EXOSAT data (Kuulkers 1995) 
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suggests that the Z sources can be divided into two 
subclasses: the Cyg-like sources (Cyg X-2, GX 5—1, 
and GX 340+0) and the Sco-like sources (Sco X-l, 
Sco X-2, and GX 17+2). Although the Z tracks 
of these sources are similar, in CDs the horizontal 
branches of the Sco-like sources are shorter and more 
inclined, while in HCDs their flaring branches extend 
toward higher countrates. Also, no HBO has been 
detected in Sco X-l or Sco X-2, while the logarithmic 
derivative of the centroid frequency with respect to 
countrate is ~ 4 for the weak HBO seen in GX 17+2, 
^3 times larger than in the Cyg-like sources (see 
Hasinger & van der Klis 1989). 

Our spectral model qualitatively reproduces the 
spectral behavior of the Sco-like Z sources if E max is 
always <0.5 keV in these sources, rather than ~0.6- 
1.0 keV as in the Cyg-like sources, since the shape of 
the spectrum in the EXOSAT energy range is then in- 
sensitive to changes in E maK . The increase in t" s f with 
increasing M on the horizontal branch then causes the 
branch to slope downward to the right in CDs, while 
the increase of T" cc and T^f on the flaring branch 
causes the countrate in the EXOSAT energy range 
to increase, creating a flaring branch that extends to 
higher countrates than in the Cyg-like Z sources. 

The lower value of E max needed to reproduce the Z 
tracks of the Sco-like sources implies that their mag- 
netic fields are smaller than the magnetic fields of the 
Cyg-like sources. This is consistent with the fact that 
no HBO has been detected in Sco X-l or Sco X-2 and 
that the HBO observed in GX 17+2 is very weak. 
Moreover, according to the unified model the log- 
arithmic derivative of the HBO centroid frequency 
with countrate should be greater for neutron stars 
with lower magnetic fields (Ghosh & Lamb 1992), in 
agreement with the behavior of the HBO frequency 
observed in GX 17+2. 

Finally, we mention that the spectral model de- 
scribed here reproduces the observed X-ray colors of 
the atoll sources if their magnetic fields are as weak 
as those of the Sco-like Z sources or weaker, and their 
mass accretion rates are ~ 0.01-0.1 Me- 

A more detailed account of our spectral model and 
its application to the Cyg-like and Sco-like Z sources 
and other accreting neutron stars in low-mass systems 
will be presented in a forthcoming paper. 
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Fig. 1. — Solid lines: Values of the spectral param- 
eters E max , y HCC , and T*f used in calculating the Z 
tracks shown in Fig. 2. Dashed lines: Values of i? max 
and y HCC computed assuming, as an example, that 
the stellar magnetic field is 7 x 10 9 G, that the radius 
of the magnetosphere as a function of M is given by 
disk-magnetosphere interaction model 2B of Ghosh 
& Lamb (1992), and that the electron density in the 
cyclotron-radiating region is constant on the horizon- 
tal branch, increases from 10 19 to 8 x 10 19 cm~ 3 down 
the normal branch (due to the rapid reduction of the 
inflow velocity caused by the increase in the outward 
radiation force), and increases a further 2% on the 
flaring branch. The bottom panel shows the accretion 
rate inferred from the spectral model and EXOSAT 
countrates, in units of the accretion rate Mi at rank 
number 2; in the physical model, Mi = Ms- 
Fig. 2. — (a) Comparison of the color-color track pro- 
duced by the spectral parameter variations shown in 
Fig. 1 with Cyg X-2 data from EXOSAT. The colors 
are those used by Hasinger & van der Klis (1989). 
(b) Similar comparison of hardncss-countrate tracks. 
As usual for EXOSAT data, the flaring branch lies 
almost on top of the normal branch in the hardness- 
countrate diagram, (c) Comparison of the predicted 
countrate spectrum at rank number 1.5 with a time- 
averaged EXOSAT countrate spectrum from near the 
middle of the normal branch. 
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